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Abstract 

We study the dynamic critical behavior of the multi-grid Monte Carlo (MGMC) 
algorithm with piecewise-constant interpolation and a W-cycle, applied to the one- 
^vq dimensional 0(4)-symmetric nonlinear cr-model [= SU(2) principal chiral model], on 

lattices from L = 128 to L = 16384. Our data for the integrated autocorrelation 
time T int ^ M 2 are well fit by a logarithmic growth. We have no idea why the critical 
1/-} slowing-down is not completely eliminated. 



By now it is widely recognized [1,2,3,4] that better simulation algorithms, with strongly 
reduced critical slowing-down, are needed for high-precision Monte Carlo studies of statistical- 
mechanical systems near critical points and of quantum field theories (such as QCD) near 
the continuum limit. One promising class of such algorithms is multi-grid Monte Carlo 
(MGMC) [5,6,7,8]. 

In the first paper of this series [6], we explained the conceptual foundations of the 
MGMC method, and proved the absence of critical slowing-down for Gaussian (free-field) 
MGMC with piecewise-constant interpolation and a W-cycle. We 1 conjectured [6, Section 
IX. C] that for asymptotically free theories — which are "almost Gaussian" except at 
very long distances — the critical slowing-down should likewise be completely eliminated 
except for a possible logarithm. However, when we made a numerical study of MGMC 
in the two-dimensional 0(iV)-symmetric nonlinear cr-model (also called iV-vector model) 
with N = 4 [8], we found to our surprise that the dynamic critical exponent is not zero: 
z int,M 2 = 0.60±0.07 (subjective 68% confidence interval). We then produced [8, Section 4.1] 
a heuristic explanation of this fact, based on the logarithmically decaying deviations from 
Gaussianness in an asymptotically free theory. This explanation does not, unfortunately, 
give a quantitative prediction for zmgmc] it merely makes the weak prediction that < 

ZMGMC < Zhe at — bath • 

Richard Brower (private communication) then pointed out to us that if our explanation 
is correct, it follows that for a one-dimensional cr-model — in which the deviations from 
Gaussianness decay much more rapidly, namely as a power law — the critical slowing-down 
should be completely eliminated, i.e. we should have zmgmc = (with no logarithm). We 
hardly doubted that this would be the case, but we decided to make an empirical test to 
settle the matter once and for all. 

In this paper we report evidence which suggests that in fact the critical slowing-down 
is not completely eliminated in the one-dimensional 0(4)-symmetric cr-model. More pre- 
cisely, we have simulated this model on periodic lattices of size L = 128 through L = 16384, 
using MGMC with piecewise-constant interpolation and a W-cycle, and we find that the 
autocorrelation time is well fit by a logarithmic growth Ti nt ^2 A + BlogL. (For L > 512 
the data can also be well fit by a power-law growth with exponent z,- nt ^2 0.18 ± 0.01, 
but the logarithmic fit is better.) Of course, it is perfectly possible that the autocorrelation 
time will level off at larger lattice sizes — even at L = 16384 it has only reached 5 sweeps 
— but the data for L < 16384 show no tendency in this direction. We are thus left with 
the distinct impression that the critical slowing-down is not completely eliminated in this 
model — even if only a logarithm remains — and we have no idea why. 

The model in question is defined by the Hamiltonian (= Euclidean action) 

L 

H((t) = — (3 (4) ^2 a x ' a x + 1 ' (1) 

x=l 

where each spin cr x is a unit vector in R 4 , and /?o(4) > is the inverse temperature in the 
0(4) normalization; periodic boundary conditions (<t_l+i = ffi) are employed. Since the 

1 Here "we" refers to the second author only. The first author is not responsible for the second author's 
sordid past. 
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sphere S is isometric (as a Riemannian manifold) to the group SU(2) via the map 



cr° + zcr 3 cr 2 + la 1 
-a 2 + 1a 1 cr° - ia 3 



(2) 



this model can also be thought of as the SU(2) principal chiral model 



H(g) = -/3 S u(2)J2 RetT (9l9x+i) 

x=l 



(3) 



with f3su(2) = /3o(4)/2. Our MGMC algorithm employs this SU(2) representation, and 
is described in detail in [8]. We used a W-cycle with one heat-bath pre-sweep and one 
heat-bath post-sweep (mi = m 2 = 1); the coarsest grid always has two sites. All runs used 
a random initial configuration ("hot start"). 

We measured static quantities (expectations) and dynamic quantities (autocorrelation 
functions) for the following observables: 



(4) 



£ 



E 



D 2Trix/L 



E 



cr ■ cr 



x+l 



(5) 
(6) 



The mean values of these observables give information on different aspects of the 2-point 
function 



G(x) = (cr -cr x ) 

G(p) = J2 eW ' X ( (T 0-(T x ) 

X 

In particular, we are interested in the susceptibility 

X = G(0) = L-'iM 2 ) 
and the analogous quantity at the smallest nonzero momentum 

F = G(p)\ M=2v/L = L-\r) . 

By combining these we can obtain the second-moment correlation length 

f( X /F)-l\ 1/2 



^4sin 2 (7r/ J L) / 



(7a) 
(7b) 



(8) 



(9) 



(10) 



The dynamic quantities for an observable A (= Ai 2 } J 7 or £) are defined from the 
unnormalized autocorrelation function 



C A (t) = (A s A s+t ) - (A)' 



(11) 
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where expectations are taken in equilibrium, and the corresponding normalized autocorre- 
lation function 

PA (t) = C A (t)/C A (0) . (12) 
We then define the integrated autocorrelation time 

2 CO 

T m t,A = ~ P^( f ) 
Z t=-oo 
-| oo 

= ? + E^W- (i3) 
z t=i 

[The factor of | is purely a matter of convention; it is inserted so that t;^^ ~ r if 
/9^(t) e~'*'/ T with r ^> 1.] We also define the exponential autocorrelation time 

T exp ,A = lim sup — - — ^ . (14) 

|*| — oo -log|^(<)| 

The integrated autocorrelation time controls the statistical error in Monte Carlo measure- 
ments of (A), while the exponential autocorrelation time measures the "slowest mode" in 
the system (provided that it is not orthogonal to A). We estimated the integrated auto- 
correlation time by standard procedures of statistical time-series analysis [9,10], using a 
self-consistent truncation window of width CTi ntj A with c = 8: see [1, Section 3] and [11, 
Appendix C]. The choice c = 8 is reasonable whenever the autocorrelation function pA(t) 
decays roughly exponentially, a behavior that we will confirm explicitly here. We estimated 
the exponential autocorrelation time by the rather subjective procedure of looking for a 
plateau in the estimates of 1/ log[/9^(t)/ p A {t + 1)] as a function of t; the error bars are also 
subjective. Of course, all estimates of T exPjA are really only approximate lower bounds on 
T~exp,A] it is impossible to rule out the existence of a (very-)slowly-decaying component of 
(very) small amplitude. 

The results of our runs are shown in Table 1. The exact finite- volume values of \ and 
£ are obtained from [12]. The deviation from the exact values has a chi-squared of 27.75 
(28 DF, confidence level = 47.8%) for the susceptibility, and 31.58 (28 DF, confidence 
level = 29.2%) for the correlation length. 2 The energy E is not shown here, but it also 
agrees with the exact values (chi-squared = 24.12, 28 DF, confidence level = 67.5%). The 
autocorrelation time of the energy, r 8 ' ntj £, is uniformly less than 0.6, so we do not bother 
to report it. 

In Figure 1 we plot Ti nt ^2 versus ^(L)/L } indicating the different lattice sizes by sym- 
bols. Clearly Ti nt ^2 is an increasing function of L at each fixed ^(L)/L } and there is no 
evidence that the increase is slowing down as L grows. 

In Figure 2 we show a log-log plot of Ti nt ^2 versus L for the data at f3 / L = 40/128 = 
0.3125 (corresponding to £(L)/L ~ 0.2). In Table 2 we report the results of a weighted 
least-squares fit to Ti nt ^2 = A L z * nt - M2 for these data, using lattice sizes L > L m i n and 
trying various values of L m i n . For L m i n = 128,256 the estimates are clearly afflicted with 
corrections to scaling: this can be seen in the chi-squared values, as well as in the clear 



2 Confidence level is the probability that the chi-squared would equal or exceed the observed value, as- 
suming that the underlying statistical model ("null hypothesis") is correct. 
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trend in the estimates of z int ^ M 2 and A; it is also apparent in the plot. For L m i n = 512 
there may be corrections to scaling on the borderline of statistical significance: the point 
at L = 512 is about 2a below the trend-line of the larger lattices, with a deviation of the 
same sign as for L = 128, 256. Being conservative, we therefore take L m i n = 1024 as our 
preferred fit: this yields the estimate z,- nt ^2 = 0.174 ± 0.011 (68% confidence interval). 
Less conservative readers might take L m i n = 512, yielding z int ^ M 2 = 0.184 ± 0.008. 

This very low (though nonzero) estimate for z int ^ M 2 suggests that the true behavior 
might be logarithmic. This suggestion is supported by the fact that in Figure 1 the 
points at each fixed £(L)/L are roughly equally spaced. In Figure 3 we show a semi- 
log plot of Ti ntt M2 versus L for the data at f3 / L = 0.3125, and in Table 3 we show the 
corresponding weighted least-squares fit to r int ^ M 2 = A + BlogL. The quality of the fit is 
much better than for the pure power-law Ansatz: even for L m i n = 128 the chi-squared is 
extremely small, and we obtain the estimates A = —1.456 ± 0.090, B = 0.6623 ± 0.0145 
(68% confidence intervals). We are left with the distinct impression that Ti nt ^2 indeed 
grows logarithmically. Of course, a definitive test of this conjecture would require lattices 
considerably larger than those used here (say, up to L ~ 10 6 ). Examination of Figure 1 
suggests that a similar logarithmic growth holds at each fixed ^(L)/L } but the constants 
A and B clearly depend on £(L)/L. 

Next we test the dynamic finite-size-scaling Ansatz 3 

r int , A (P,L) « Z(P,Ly^g A (ttP,L)/L) (15) 

for A = Ai 2 . Here g A is an unknown scaling function, and g A (0) = fim^o 9a{x) is supposed 
to be finite and nonzero. 4 We impose z int ^ M 2 = 0.174 — which is a reasonable "effective 
exponent" for our range of L — and plot r int } M 2 C~ Z ' nt ' M2 versus ^(L)/L; the points should 
fall on a single curve, at least asymptotically for large L. The result is shown in Figure 4. 
There are small corrections to scaling for L < 512, just as observed in Table 2; but the 
points do clearly appear to be approaching a limiting curve. 

Finally, we want to test the more detailed dynamic finite-size-scaling Ansatz 

p A (t;/3,L) w \t\-^h A (t/r expA (j3,L); Z(/3,L)/L) , (16) 

3 Our preceding results suggest that we ought to use instead a logarithmic Ansatz, e.g. 

n n t,A(P>L) « \[og£{p, L) + c{£{p,L) I L)]g A {t{p,L) I L) . 

But it is difficult to know what to use for the function c(£(/?, L)/L). If one simply uses c = 0, the agreement is 
poor; clearly the additive constant cannot be neglected for our range of L, as is obvious already from Figures 1 
and 3 and Table 3. So we decided to use instead the power-law Ansatz as a reasonable "phenomenological" 
fit. 

4 It is of course equivalent to use the Ansatz 

TimAP^) ~ L z ^g A (c(j3,L)/L^ , 

and indeed the two Ansatze are related by 'cia(x) = x z ' nt - A gA(x). However, to determine whether 
lim^jo 9a(%) = linicjo x~ z ' nt - A 'gA(%) is nonzero, it is more convenient to inspect a graph of qa than one 
of §a- 
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where pa is an unknown exponent and Iia is an unknown scaling function. Summing (16) 
over t, it follows that 5 

^m,A ~ r\~H . (17) 
Therefore (16) can equivalently be rewritten as 

p A (t;P,L) « \t\- pA h A [tlr int , A {^Lfl^-^; Z(P,L)/L) . (18) 

In Figure 5(a) we test this latter Ansatz for A = Ai 2 } using the inspired guess pj^2 = 0, 
for the data at f3 / L = 0.3125. The data fall beautifully on a single curve (until the signal 
is obscured by the statistical noise), in excellent agreement with (18). 

Thus, although pa = does not hold in general, it does appear empirically to hold 
(at least for the global observable A = Ai 2 ) in our one-dimensional model. In such a 
situation, the ratio Ti nt ^2 /r exp ^2 tends as L — > oo to a limiting value (which could depend 
on i{L)l ' L) that is characteristic of the dynamic universality class. We can estimate this 
value by a linear fit 

\ogp M 2(t) = a-bt/T mtM 2 (19) 

to the data in Figure 5(a), restricting attention to the interval 1.5 < t/Ti nt ^2 < 3 where the 
plot indeed appears linear; the slope b gives the ratio Ti nt ^2 /r exp ^2. We obtain a = —0.18, 
b = 0.90. Of course, as noted above, all estimates of T exp ^2 are really only approximate 
lower bounds; so the true ratio r int ^ M 2 / r exp ^ M 2 might be lower than this. 

In Figure 5(b) we make an identical plot, but now using the data from all values of 
f3 / L > 0.3125. (The data at smaller f3 / L are omitted because their statistics are much 
poorer.) Surprisingly, the points continue to fall well on a single curve; the linear fit (19) 
on the interval 1.5 < t/Ti nt ^2 < 3 yields a = —0.16, b = 0.92. This indicates that the 
dependence on £(L)/L in (16) / (18), if any, is extremely weak. 

We also attempted to estimate r exp ^ M 2 by looking for a plateau in the estimates of 
1 1 \og[pM. 2 {t) I PM 2 {t + 1)] as a function of t. More specifically, we used the largest t before 
the values of 1/ log[pj^2(t)/ pj^2(t + 1)] begin to decrease or oscillate significantly. This 
roughly corresponds to choosing t by any one of the following three criteria 6 : 

(a) t l.lT mtM 2 

(b) t such that PM 2 (t) ~ 0.17 

(c) t such that Pm 2 {^) ~ 50 [^-^^2 /(Sweeps — Discard)] 1 / 2 

The error bars are subjective 68% confidence intervals. The resulting estimates of T exp ^ M 2 } 
for the data at f3 / L = 0.3125, are given in Table 4. We see that the ratio r int ^ M 2 / r exp ^ M 2 
is essentially independent of L and equals 0.90. This agrees with the estimate from 
Figure 5(a). Again, all estimates of T exp ^2 are only approximate lower bounds. 

5 Note that, contrary to much belief, Ti ntt A and T exPt A do not in general have the same dynamic critical 
exponent z; rather, we have from (17) the dynamic scaling law z inttA = (I— Pa) Zexp,A- Thus, Z{ n -t^A — ^exp,A 
only if pa = 0. For further discussion, see [1,4]. 

6 All three criteria translate, in different ways, the intuitive idea of using the largest t such that the signal- 
to-noise ratio is "not too small". In particular, criterion (c) comes from the fact that the statistical error 
bar on Pa{^) is roughly of order [n n t t A/ (Sweeps — Discar d)] 1 / 2 : see [1,9,10,11]. 
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d = 1 0(4) model 


T 

±j 


P 


^ 117/^/^ C 1 1 1 1 C ^ 1* / \ 

O Wccpo JJloCdl U. 


X Xexact 


S> Xexact 


T int,M 2 


128 


10 


50000 5000 

\j \j \j \j \j \j \j 


12.74 ( 0.04) 12.72 


6 41 ( 15) 6 34 


59 (0 01) 


128 


20 


42869 5000 


26.03 ( 0.10) 26.02 


13.05 ( 0.13) 13.00 


0.93 (0.03) 


128 


40 


200000 5000 


51.41 ( 0.08) 51.41 


26.13 ( 0.07) 26.14 


1.75 (0.03) 


128 


60 


200000 5000 


70.51 ( 0.08) 70.56 


37.94 ( 0.08) 37.96 


1.97 (0.04) 


128 


80 


200000 5000 


82.93 ( 0.06) 82.89 


47.84 ( 0.09) 47.78 


1.81 (0.03) 


128 


100 


200000 5000 


90.99 ( 0.05) 90.98 


56.02 ( 0.09) 56.02 


1.65 (0.03) 


128 


120 


200000 5000 


96.64 ( 0.04) 96.63 


63.24 ( 0.09) 63.21 


1.58 (0.03) 


256 


20 


50000 5000 


25 86 < 09) 26 02 


1 2 42 < 32) 1 3 00 


70 (0 02) 


256 


40 


50000 5000 


52.60 ( 0.22) 52.67 


26.28 ( 0.27) 26.33 


1.23 (0.04) 


256 


80 


200000 5000 


103.30 ( 0.19) 103.38 


52.56 ( 0.16) 52.60 


2.21 (0.04) 


256 


120 


200000 5000 


141.54 ( 0.17) 141.48 


76.26 ( 0.18) 76.18 


2.31 (0.04) 


256 


160 


200000 5000 


166.04 ( 0.14) 166.02 


95.84 ( 0.19) 95.78 


2.17 (0.04) 


256 


200 


200000 5000 


182.22 ( 0.11) 182.12 


112.39 ( 0.20) 112.24 


2.06 (0.04) 


256 


240 


200000 5000 


193.17 ( 0.10) 193.37 


126.11 ( 0.21) 126.58 


1.98 (0.04) 


512 


40 


31080 5000 


52.78 ( 0.27) 52.68 


27.07 ( 0.84) 26.33 


0.82 (0.03) 


512 


80 


31159 5000 


106.51 ( 0.61) 105.98 


53.85 ( 0.75) 53.00 


1.42 (0.06) 


512 


160 


400000 5000 


207.57 ( 0.29) 207.32 


105.66 ( 0.25) 105.51 


2.71 (0.04) 


512 


240 


200000 5000 


283.98 ( 0.37) 283.34 


153.25 ( 0.39) 152.64 


2.77 (0.06) 


512 


320 


200000 5000 


331.65 ( 0.31) 332.27 


190.97 ( 0.41) 191.78 


2.62 (0.05) 


512 


400 


200000 5000 


364.43 ( 0.25) 364.41 


224.58 ( 0.43) 224.67 


2.44 (0.05) 


512 


480 


200000 5000 


386.93 ( 0.21) 386.86 


253.53 ( 0.45) 253.33 


2.32 (0.05) 


1024 


80 


44399 5000 


106.00 ( 0.46) 106.00 


55.28 ( 1.39) 53.00 


0.91 (0.02) 


1024 


160 


50000 5000 


213.40 ( 1.01) 212.62 


106.98 ( 1.25) 106.33 


1.67 (0.06) 


1024 


320 


400000 5000 


414.81 ( 0.62) 415.21 


211.02 ( 0.54) 211.34 


3.13 (0.05) 


2048 


640 


400000 5000 


832.50 ( 1.33) 830.99 


424.63 ( 1.15) 423.00 


3.58 (0.06) 


4096 


1280 


400000 5000 


1665.37 ( 2.83) 1662.55 


849.15 ( 2.46) 846.31 


4.05 (0.07) 


8192 


2560 


200000 5000 


3308.68 ( 8.45) 3325.66 


1678.94 ( 7.34) 1692.94 


4.43 (0.12) 


16384 


5120 


100000 5000 


6614.63 (25.93) 6651.90 


3354.28 (22.46) 3386.19 


5.09 (0.21) 



Table 1: The results of our runs for the one-dimensional 0(4) model. "Sweeps" is the total 
number of MGMC iterations performed; "Discard" is the number of iterations discarded 
prior to beginning the analysis. Error bar (one standard deviation) is shown in parentheses. 
The exact values are obtained from [12]. 
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d = 1 0(4) 


model 








z int ,A4 2 


A 


x 2 


128 


0.223 ( 0.005) 


0.6396 ( 0.0219) 


53.264 ( 


6 DF, level = 


0.000 %) 


256 


0.200 ( 0.006) 


0.7647 ( 0.0344) 


15.758 ( 


5 DF, level = 


0.757 %) 


512 


0.184 ( 0.008) 


0.8715 ( 0.0504) 


2.834 ( 


4 DF, level = 


58.602 %) 


1024 


0.174 (0.011) 


0.9460 (0.0829) 


1.285 (3 DF, level = 


73.276 %) 


2048 


0.163 ( 0.017) 


1.0352 ( 0.1459) 


0.618 ( 


2 DF, level = 


73.431 %) 


4096 


0.155 ( 0.029) 


1.1144 ( 0.2787) 


0.490 ( 


1 DF, level = 


48.378 %) 


8192 


0.200 ( 0.071) 


0.7283 ( 0.4784) 


0.000 ( 


DF, level = 


100.000 %) 



Table 2: Weighted least-squares fit for r mt)M 2 = AL z ' nt - M2 at /3/L = 40/128 = 0.3125, 
using lattice sizes L > L m i n . Errors in parentheses represent one standard deviation, and 
"DF" stands for "degrees of freedom". Confidence level is the probability that % 2 would 
equal or exceed the observed value, assuming that the underlying statistical model ("null 
hypothesis") is correct. The line in boldface marks our preferred fit. 
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d = 1 0(4) 


model 








A 


B 


x 2 


128 


-1.456 (0.090) 


0.6623 (0.0145) 


1.678 (6 DF, level = 


94.681 %) 


256 


-1.410 ( 0.139) 


0.6560 ( 0.0207) 


1.493 ( 


5 DF, level = 


91.392 %) 


512 


-1.310 ( 0.201) 


0.6428 ( 0.0282) 


1.016 ( 


4 DF, level = 


90.737 %) 


1024 


-1.446 ( 0.336) 


0.6598 ( 0.0438) 


0.759 ( 


3 DF, level = 


85.931 %) 


2048 


-1.512 ( 0.593) 


0.6676 ( 0.0729) 


0.741 ( 


2 DF, level = 


69.042 %) 


4096 


-1.585 ( 1.153) 


0.6759 ( 0.1341) 


0.735 ( 


1 DF, level = 


39.113 %) 


8192 


-4.150 ( 3.205) 


0.9522 ( 0.3489) 


0.000 ( 


DF, level = 


100.000 %) 



Table 3: Weighted least-squares fit for r mt)M 2 = A + B log L at /3/L = 40/128 = 0.3125, 
using lattice sizes L > L m i n . Errors in parentheses represent one standard deviation, and 
"DF" stands for "degrees of freedom". Confidence level is the probability that % 2 would 
equal or exceed the observed value, assuming that the underlying statistical model ("null 
hypothesis") is correct. The line in boldface marks our preferred fit. 
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d = 1 0(4) mo 


del 


L 


T int,M 2 




T int ,M 2 / T exp ,M 2 


128 


1.75 (0.03) 


1.9 (0.2) 


0.92 


256 


2.21 (0.04) 


2.4 (0.3) 


0.92 


512 


2.71 (0.04) 


3.1 (0.2) 


0.87 


1024 


3.13 (0.05) 


3.4 (0.2) 


0.92 


2048 


3.58 (0.06) 


4.0 (0.3) 


0.90 


4096 


4.05 (0.07) 


4.5 (0.3) 


0.90 


8192 


4.43 (0.12) 


4.9 (0.4) 


0.90 


16384 


5.09 (0.21) 


5.8 (0.4) 


0.88 



Table 4: Estimates of Ti nt ^2 and T exv ^ 2 for the points with f3 / L = 0.3125. See text 
the method of calculation. Errors in parentheses represent one standard deviation. 
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d=l 0(4) model, MGMC W-cycle 



4 



2 







0.0 



CD 

□ 



I 
1 



0.2 



£(L)/L 



0.4 



m 



Figure 1: Plot of Ti nt ^2 versus £(L)/L for the W-cycle MGMC algorithm, for lattice sizes 
L = 128 (□), L = 256 (O), L = 512 (O), L = 1024 (*), L = 2048 (+), L = 4096 (x), 
L = 8192 (hfn), L = 16384 (&<)• 
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T int m 2 vs - L: Power — law fit 
10 i 1 — i i i i m i 1 — r— 




10 2 10 3 10 4 10 5 

L 



Figure 2: Test of power-law Ansatz r int ^ M 2 = A L z * nt - M2 } for data at f3 / L = 0.3125. The 
straight line is the least-squares fit T int ^ M 2 = 0.9460 L 0A74 , from Table 2 with L m i n = 1024. 
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Fi gure 3: Test of logarithmic Ansatz Tj^j ^2 = A + BlogL, for data at f3/L = 0.3125. The 
straight line is the least-squares fit Ti nt ^2 = —1.456 + 0.6623 logL, from Table 3 with 

Lfnin — 128. 
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Figure 4: Finite-size-scaling plot of r int ^/f'"'^ 2 versus £(L)/L for the W-cycle MGMC 
algorithm, for lattice sizes L = 128 (□), L = 256 (O), L = 512 (O), £ = 1024 (*), 
L = 2048 (+), L = 4096 (x), L = 8192 (*), L = 16384 (&). Here z intjM 2 = 0.174. 
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Figure 5: Autocorrelation function p M 2 as a function of time rescaled by T int ^ M 2 } for the 
case of (a) f3 / L = 0.3125, and (b) all points with f3 / L > 0.3125. The straight lines are 
linear fits to the interval 1.5 < t/ri nt ^2 < 3. 



16 



